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I. Basis of the report 

1 . With regard to the elements of the international application (Replacement sheets which have been furnished to 
the receiving Office in response to an invitation under Article 14 are referred to in this report as "originally filed" 
and are not annexed to this report since they do not contain amendments (Rules 70. 16 and 70. 17)): 
Description, pages: 

2,3 as originally filed 

1 ,4-1 0 as received on 20/08/2001 with letter of 20/08/2001 
Claims, No.: 

1-10 as received on 20/08/2001 with letter of 20/08/2001 
Drawings, sheets: 

1/7-7/7 as originally filed 

2. With regard to the language, all the elements marked above were available or furnished to this Authority in the 
language in which the international application was filed, unless otherwise indicated under this item. 

These elements were available or furnished to this Authority in the following language: , which is: 

□ the language of a translation furnished for the purposes of the international search (under Rule 23.1 (b)). 

□ the language of publication of the international application (under Rule 48.3(b)). 

□ the language of a translation furnished for the purposes of international preliminary examination (under Rule 
55.2 and/or 55.3). 

3. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the 
international preliminary examination was carried out on the basis of the sequence listing: 

□ contained in the international application in written form. 

□ filed together with the international application in computer readable form. 

□ furnished subsequently to this Authority in written form. 

□ furnished subsequently to this Authority in computer readable form. 

□ The statement that the subsequently furnished written sequence listing does not go beyond the disclosure in 
the international application as filed has been furnished. 

□ The statement that the information recorded in computer readable form is identical to the written sequence 
listing has been furnished. 

4. The amendments have resulted in the cancellation of: 
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□ the description, pages: 

□ the claims, Nos.: 

□ the drawings, sheets: 



5. □ This report has been established as if (some of) the amendments had not been made, since they have been 
considered to go beyond the disclosure as filed (Rule 70.2(c)): 

(Any replacement sheet containing such amendments must be referred to under item 1 and annexed to this 
report.) 



6. Additional observations, if necessary: 



V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

1. Statement 

Novelty (N) Yes: Claims 1-10 

No: Claims 

Inventive step (IS) Yes: Claims 1-10 

No: Claims 

Industrial applicability (I A) Yes: Claims 1-10 

No: Claims 



2. Citations and explanations 
see separate sheet 



VII. Certain defects in the international application 

The following defects in the form or contents of the international application have been noted: 
see separate sheet 
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Re Item V 

Reasoned statement under Article 35(2) with regard to novelty, inventive step or 
industrial applicability; citations and explanations supporting such statement 

1 . Reference is made to the following document: 

D1: J.D. WALLIS et al.: "An Optimal Rotator for Iterative Reconstruction", IEEE 
TRANSACTIONS ON MEDICAL IMAGING, IEEE INC., US, NEW YORK, vol. 
16, no. 1, 1 February 1997, pages 118-123, XP000685496 ISSN: 0278-0062 

2. The preferred use of the techniques described in D1 is image rotation (inverse or 
forward mapping are discussed). In contrast, the main purpose described in the 
present application is image enlargement. However, the claimed subject-matter is not 
limited to enlargement but has wider applications, and is in fact applicable whenever 
in-between target values must be interpolated from a one- or multi-dimensional grid 
of discrete source values (cf. page 1 , lines 9-1 1 , of the application). 

Therefore, despite of their different contexts, D1 discloses relevant prior art. 
However, D1 (page 118, abstract, last 2 lines of the left column, first 14 and last 2 
lines of the right column, page 119, first partial paragraph) only represents 
background prior art as it discloses performing Gaussian weighted interpolation 
amongst initial pixel values so as to determine pixel values at fractional pixel 
locations. 

Starting from these known techniques, the present invention allows to adjust the 
value obtained by (Gaussian) weighted interpolation such that edges are preserved 
but excessive blurring is avoided. 

This is allowed, in particular, by the claimed determination of a measure of dynamics 
in a local region around the target (pixel) value to be interpolated. 
Neither this claimed measure of dynamics, nor the final claimed adjustment are 
disclosed or even only suggested in D1 or in the other documents cited in the 
international search report. 

Consequently, notwithstanding the objection of Section VIII below, the subject-matter 
set out in the present claims, and particularly in claim 1 , is considered to be novel 
and non-obvious with respect to the disclosures of the available prior art. It is also 
evident that the invention is industrially applicable. 

The requirements of paragraphs (1) to (4) of Article 33 PCT are thus met. 
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Re Item VII 

Certain defects in the international application 

3. In claim 1 , lines 6-7, "the target value is calculated by weighted interpolation is 
adjusted" is linguistically inconsistent and should instead read "the target value as 
calculated by weighted interpolation is adjusted". 

Correspondingly, at line 3, "interpolation" should read "weighted interpolation", to 
create an antecedent for the "as calculated by weighted interpolation" of lines 6-7. 
In claim 1 the expression "measure of the dynamics" should be replaced by "measure 
of dynamics" for linguistic reasons. For consistency, also the dependent claim should 
use this expression instead of the variety of similar expressions used. 
In dependent claim 7, " as second region" should read "the local region", and "which 
second region is optionally of the same size as the first region" should read "which 
local region coincides with the local region" for consistency with claim 1 and with the 
description, cf. page 7, lines 16-18 (page 7, lines 8-10 of the originally filed 
description). 

Dependent claim 8 appears superfluous with respect to dependent claim 2. 
There is no antecedent in the preceding claims for "the average" mentioned at claim 
9, line 1. 

4. The opening part of the description (cf. bottom of page 2) should have been modified 
to bring it into agreement with the amended claims, Rule 5.1 (a) (iii) PCT. 

In order to meet the requirements of Rule 5.1 (a) (ii) PCT, document D1 should have 
been acknowledged and briefly discussed in the opening part of the description. 
In the description, the same text is repeated at lines 25-31 of page 3 and lines 1-7 
(till "the second region") of page 4. 
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International application No. 

PCT/NL 00/00359 


International filing date (day/month/year) 

25/05/2000 


(Earliest) Priority Date (day/month/year) 

31/05/1999 


Applicant 

ESSENTIAL RESOURCE B.V. 



This International Search Report has been prepared by this International Searching Authority and is transmitted to the applicant 
according to Article 1 8. A copy is being transmitted to the International Bureau. 

This International Search Report consists of a total of 3 sheets. 

f"X~| It is also accompanied by a copy of each prior art document cited in this report. 



1 . Basis of the report 

a. With regard to the language, the international search was carried out on the basis of the international application in the 
language in which it was filed, unless otherwise indicated under this item. 

I I the international search was carried out on the basis of a translation of the international application furnished to this 
Authority (Rule 23.1 (b)). 

b. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the international search 
was carried out on the basis of the sequence listing : 
| | contained in the international application in written form. 

filed together with the international application in computer readable form, 
furnished subsequently to this Authority in written form, 
furnished subsequently to this Authority in computer readble form. 



2. 
3. 



□ 
□ 
□ 
□ 

□ 

□ 
□ 



the statement that the subsequently furnished written sequence listing does not go beyond the disclosure in the 
international application as filed has been furnished. 

the statement that the information recorded in computer readable form is identical to the written sequence listing has been 
furnished 

Certain claims were found unsearchable (See Box I). 
Unity of Invention Is lacking (see Box II). 



4. With regard to the title, 

[X] the text is approved as submitted by the applicant. 

| | the text has been established by this Authority to read as follows: 



5. With regard to the abstract, 

|X| the text is approved as submitted by the applicant. 

□ the text has been established, according to Rule 38.2(b), by this Authority as it appears in Box III. The applicant may, 
within one month from the date of mailing of this international search report, submit comments to this Authority. 

6. The figure of the drawings to be published with the abstract is Figure No. 3 

[X] as suggested by the applicant. None of the figures. 

| | because the applicant failed to suggest a figure. 

| | because this figure better characterizes the invention. 
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Category " Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



US 5 561 724. A 
1 October 1996 
abstract 
column 6, line 
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25 - line 35 
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WALLIS J W ET AL: "AN OPTIMAL ROTATOR FOR 
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IMAGING, US, IEEE INC. NEW YORK, 
vol . 16, no. 1, 
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"O" document referring to an oral disclosure, use, exhibition or 
other means 

*P* document published prior to the international filing date but 
later than the priority date claimed 



T" later document published after the intemationaJ filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y* document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 
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abstract 
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abstract 
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Method for data processing 

The present invention relates to a method for processing data in the form of a grid of discrete 
source values, wherein at least one target value within a region of source values is deteriuined 
by means of interpolation in said region of source values. 

Such a method can be used for data of diverse nature, and particularly to expand, compress 
or decompress computer data with a non-critical bit value, such as for instance sound and 
animation files. The starting point is always a one- or multi-dimensional grid of discrete source 
values among which concrete target values have to be predicted by means of interpolation. 
The source values herein comprise numerical values which, in the case of for instance an 
image, indicate the colour intensity of a basic colour present therein and, in the case of a 
sound file, represents the frequency, i.e. the pitch, or the intensity thereof. In the case of other 
types of data a different information component in the source data will similarly be taken as 
measure. 

For digital images there exists a number of more or less standardized forms in which the data 
is stored. The so-called bitmap format is for instance a form wherein the data is stored in a 
series of discrete numerical values. This series can be converted comparatively simply into a 
two-dimensional matrix of source values which each indicate the colour intensity of pixels 
corresponding therewith. In a monochromatic image there thus results a single matrix, while in 
a colour image such a matrix can be constructed for each of the basic colours. 

One of the simplest ways of enlarging such an image is to expand the discrete pixels. The 
result hereof is shown in figures 1 and 4, wherein the lower image T 1 is roughly a five-fold 
enlargement of the source image S. It can be seen clearly that such a method of enlargement 
detracts significantly from the resolution of the final image. Because each pixel of the source 
image S reappears in the result Tl as a rectangular image area of a uniform colour (intensity), 
a clear block structure is created which does not correspond with the original source image. 



AMENDED SHEET 



ussioei.wo.ob2 



- — WW all* WW WW 



uwixvWjUUACiiU L1UL 



14)003 



-4- 



ln the context of the present invention local dynamics in the source values are understood to 
mean the degree in which the source values in the relevant region differ from each other and 
also the steepness with which these differences are present. A further particular embodiment 
of the method according to the invention comprises in this respect an algorithm which assigns 
a numerical value hereto and is characterized for this purpose in that the dynamics are derived 
from a normalized difference between a source value and an average of all source values in 
the second region. For the average of all source values in the second region a weighted 
average is herein preferably taken which assigns a heavier weighting to source values located 
closer in the grid than to source values located further away and which particularly utilizes 
non-linear density distribution for the purpose of determining the weighting factors and 
particularly a Gaussian distribution, at least an exponential density distribution. 



a 

more 



;in a 



In order to determine the different factors which play a part in predicting the target value 
within the scope of the invention, the starting point is always a number of source values h 
finite first and second region as frame of reference. Relatively far-removed source values can 
be ignored because they make no or hardly any contribution towards the precision with which 
the target value is calculated and can even detract from this precision in that there is alack of 
any relationship with the target value for predicting. The calculation time of the method 
moreover increases when more source values are taken as frame of reference. Conversely, 
only one or two reference values will in many cases be too few to enable prediction of the 
target value in reliable manner. A particular embodiment of the method according to the 
invention has been found effective in practice wherein the first and second region both extend 
over nine source values in the grid. 

The adjustment of the interpolated value to the local maximum or minimum within the scope of 
the present invention can be performed per se in diverse ways. A further particular 
embodiment of the method according to the invention is however characterized in that the final 
target value is a weighted average of the target value determined on the basis of interpolation 
and the local maximum and minimum, wherein a weighting factor is employed which depends 
on average local dynamics of the source values located around the target value and the 
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relative location of the target value determined on the basis of interpolation relatives the local 
maximum and minimum. 

In order to determine the local minima, maxima and dynamics it isper se possible to carry out 
a separate calculation for each target value to be determined. A preferred embodiment of the 
method according to the invention is however characterized in that for all source values an 
associated local rninimum and maximum and an associated dynamic value is determined 
beforehand so as to be read for the purpose of deterrnining the target value. A calculation is 
herein performed once for all source values together in order to determine said values, so that 
these are men immediately available. This saves considerable calculation time since 
calculations which would otherwise be performed repeatedly now take up calculation time 
once-only. 

The invention is not only suitable for data with only a single information component but also 
for data in the form of source values with separate numerical values for separate information 
components. In accordance with a further embodiment of the method according to the 
invention, for each information component a target value located among the sourcevalues is 
herein determined individually for this information component 

The invention will now be further elucidated with reference to an embodiment and an 
associated drawing. In the drawing: 

figurel shows the enlargement of a first image based on an existing technique making using of 

pixel expansion; 

%ure2 shows an enlargement of the first image based on a second existing technique making 

using solely of pixel interpolation; 
fi gure 3 shows an enlargement of the first image based on an embodiment of the 

method according to the present invention; 
figure 4 shows the enlargement of a second image based on an existing technique 
making using of pixel expansion; 
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figure 5 shows an enlargement of the second image based on a second existing 

technique making using solely of pixel interpolation; and 

figure 6 shows an enlargement of the second image based on the embodiment of the 

method according to the present invention used in figure 3. 

Although the invention lends itself in principle to any type of data wherein the exact bit value is 
not critical, it is best elucidated within the scope of the present invention on the basis of an 
example with graphic data in the form of an image S. The data on the basis of which the image 
can be constructed is possibly already in the form of a grid or matrix corresponding with the 
pixels of the image, with for each pixel a numerical source value indicating a colour value or, 
with a view to the method according to the invention, is brought into such a form. In the case 
of a monochromatic image the colour value amounts only to an intensity, for instance a grey 
value; in the case of a colour image this will be a set of for instance three values, one for each 
basic colour. A much used format in which images are digitally stored is the so-called bitmap 
format which, for each pixel of a colour image, comprises three channels of 8 bits each, and 
therefore has 256 values for each basic colour. This format can be employed relatively 
directly by the method according to the invention, other formats possibly require a conversion 
similar to that with which such a format is transferred to the (image) memory of a computer. 

Assuming a black/white image S, a part of the thus assembled matrix is shown by the grid of 
figure 7, wherein the points in the grid represent the pixels of the source image S with for each 
grid point S fJ a discrete value of ^ from 0 to 255 from the matrix to thereby designate the grey 
value of the relevant pixel. In order to enlarge such an image while retaining resolution, 
intermediate points have to be created, of which the value of the grey value has to be 
calculated. One of these points T is shown in the figure. 

A first estimate of the grey value \ of the intermediate pixel T can be obtained by interpolation 
on the basis of the values of the surrounding pixels. The point of departure herefor in this 
example is a weighted average of the values of the surrounding pixels, wherein the weighting 
assigned to each pixel is highly dependent on its distance from T. A more particular point of 
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departure herein for the weighting factors is a normalized Gaussian distribution, of which pixel 
T forms the origin. Because of the sharp fall in the weighting factors according to this curve, a 
finite region A of for instance 4x4 pixels around T suffices for the interpolation, which pixels 
are hatched in the figure. On the basis of this interpolation an interpolated grey value P t results 
5 for the pixel T to be calculated in accordance with the following relation: 

p * = 2 GjfJp, where ]T Go- 1 

According to the invention a local maximum and a local rninimum around the pixel T to be 
calculated is also determined in addition to the interpolated value P r In the embodiment 
described here the starting point for this purpose is a reference region B of 4x4 pixels S., . 
□ 1 0 i£oa~.S2j, an origin S Of0 of which is formed by the pixel in the original image lying just in front 
2 of the pixel for calculating. This region is shown hatched in the figure. A local maximum I n 

^ and a local minimum ^ of the grey values of the original pixels are determined within this 

i Q region. 



In addition, the local dynamics in the grey values are determined according to the invention 
within a reference region around the pixel T for calculating. The starting point therefor in this 
embodiment is the same reference region B as that in which the local minimum and maximum 
have been determined. These dynamics provide a measure for the hardness or contrast of the 
image and are derived in this embodiment from the normalized average difference between the 
grey values of the pixels Sy in reference region B and a weighted average of the pixels in a 
reference region C of 5x5 pixels around the pixel Sg for which the local dynamics has to be 
calculated. This region C therefore differs for each pixel S„ and is indicated by way of 
example in the figure for S,». This is expressed in the following formula: 

where u~un = £ with ~£ G v = 1 and Y G„ = 1 
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Thestarting point for the weighted average 1^ here are weighting factors G' M on the basis of 
a noimalized Gaussian distribution which attributes a greater weighting to pixels close to the 
pixel T for calculating than to further removed pixels. For the resulting factor D fj there applies: 
0 < Dy < 1 . A value D™ 1 indicates that the source image in the relevant region is very hard, 
i.e. has very sharp transitions in intensity, while a value Dij=0 indicates precisely the opposite, 
Le. no differences in intensity occur within the region. The local dynamics D t associated with 
the pixel T to be calculated are derived from the weighted average of the values D 0 - from the 
reference region B around the pixel T to be calculated. 

The grey value of the pixel T to be calculated is established by adjusting the interpolated value 
P, to the local maximum or minimum on the basis of the thus found hardness value D s . Taken 
for this purpose in this embodiment is a weighted average of P r and the local maximum and 
minimu m, wherein a weighting factor is applied which depends on the difference from P t and is 
proportional to the hardness value D t . This is expressed for instance by the following formula 
applied in this 
example: 



Im/rt + (lm<zx — Im irr) . SHI rt( 1 

\ Im ax — Im in J 



where sino(f) = J4.sim>r(r - K) and sin„(0 = sino.sin„- i(f) with n* 5 

If in the relevant reference region B the source image is completely level and has no 
differences in intensity, which is reflected by a value D=0, no shifting or adjustment of the 
interpolated value P te takes place. In the case of a comparatively hard image with great 
differences in intensity, wherein D will move in the direction of a value D= 1 3 the adjustment 
will in contrast be maximal. The final value I t will in that case be shifted to the local minimum 
or maximum, depending on which value lies closest to the interpolated value P^ thereby 
resulting in more definition. The local dynamics of the original image S are thus taken into 
account in the final grey value I D which is calculated for the pixel T to be added. 
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The foregoing algorithm is applied for all pixels which have to be calculated among or 
adjacently of the original pixels in order to create the desired enlargement while retaining 
resolution. In practice this means mat for all pixels of the original source image the associated 
local maxima, minima and average hardness values are calculated beforehand so that these are 
then immediately available for the above stated calculations. The stated algorithm allows of 
simple translation into computer software with which a suitable computer can be loaded to 
perform the calculations fully automatically. In a user interface adapted thereto setting options 
can be offered for fme-tuning of the size of the reference regions, the relations for the 
weighting factors and other parameters. 

The result on the basis of this embodiment of the method according to the invention is shown 
in figures 3 and 6 for respectively a relatively hard source image in the form of a sharply 
defined character and a more variegated source image corresponding with a photographic 
image. It can clearly be seen that the definition of the character in figure 3 has been retained 
despite the great enlargement thereof while the variegation in the image of figure 6 is retained 
with the method according to the invention and is even better preserved than on the basis of 
solely an interpolation method as shown in figure 5. 

Although the invention has been further elucidated above with reference to only two 
embodiments, it will be apparent that the invention is by no means limited to the given 
embodiments. The invention can thus also be applied to images with more colours, wherein 
the above stated algorithm is applied separately to each basic colour. The invention can 
advantageously also be applied for the interpolation of sound data and other data with 
noncritical bit values while retaining dynamics. The above stated weighting factors and frames 
of reference, although very effective, have been given solely as examples. Depending on the 
concrete situation, a fine-tuning thereof can take place to further enhance the quality of the 
final result. Larger, smaller or identical reference regions can for instance thus be applied to 
calculate the various above stated factors and, in order to determine the final grey value of the 
pixel for calculating, a different algorithm can also be chosen which takes account of the local 
dynamics derived in this pixel and the local maximum and minimum 
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The invention generally provides a method of interpolating data while retaining dynamics, 
which implies that the character and the definition of the source data is retained to a significant 
degree in the result. 
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Claims 

1 . Method for processing data in the form of a grid of discrete source values, wherein at 
least one target value (T) within a region (A) of source values is determined by means of 
interpolation in said region (A) of source values characterised in that a rninimum value (1^ 
and a maximum value are determined within a local region of source values around the 
target value (T), in mat a measure of the dynamics is determined within a local region (B) of 
source values around the target value (T) and in that the target value is calculated by weighted 
interpolation is adjusted in the direction of either said minimum value or said maximum value 
on the basis of said determined measure of the dynamics. 

2. Method according to claim 1 characterised in that said measure of dynamics is 
determined as a normalized weighted value of the absolute differences - I gcm ||) in source 
values within said local region (B). 

3. Method according to claim 2 characterized in that one of said absolute differences 
(Iffij - I g cm(ij)ll) is calculated for each one of the source values in said local region (B), and in 
that each said difference Qfly - I^H) is calculated between a given one (l ijs ijeB) of said 
source values and a weighted average (I gcm(ij) ) of source values in a further local region (C) 
corresponding to said given one source value. 

4. Method according to any of the preceding claims charactered in that the direction in 
which said adjustement is performed depends on the relanVe^erence between said target 
value calculated by weighted interpolation (PJ and saioVmlnimum and niaximum value 

5. Method according to any of the preceding claims characterized in that use is made of 
weighted interpolation on the basis of^non-linear density distribution which assigns a heavier 
weighting to source values locateo>c1oser in the grid than to source values located further 
away, in particular a Gaussian/istribution, at least an exponential density distribution. 
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which Hes in the grid closest to the targe^y^^deternnned, is taken as source of a 
region extending over a finite numbe^of mutually adjacent source values and that the local 
maximum and the local minimurrfare determined in this region. 

7. Method as claimed in claim 6, characterized in that the measure for the dynamics of 
the source values is determined in a second region extending over a finite number of mutually 
adjacent source values, which second region is optionally of the same size as the first region in 
which the local maximum and min-jr^m ^ determined. 

8. Method as claimed in claim 7, characterized in that the dynamics are derived from a 
normalized difference between a source value and an average of all source values in the 
second region. 

9. Method as claimed in claim 8, characterized in that for the average of all source values 
in the second region a weighted average is taken which assigns a heavier weighting to 

source values located closer in the grid than to source values located further away and which 
particularly utilizes a non-linear density distribution for the purpose of deterrnining the 
weighting factors and more particularly from a Gaussian distribution, at least from an 
exponential density distribution. 

10. Method according to any of thep^cetf^ in that the final target 
value is a weighted average of the tar^alu^te^rmned on the basis of interpolation and the 
local maximum and nunimum, whereina^We^ting factor is employed which depends on 
average local dynamics of the sotpee values located around the target value to be detennined 
and the relative location c^etoget value determined on the basis of interpolation relative to 
the local maximum and minimum. 
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Method for data processing 

The present invention relates to a method for processing data in the form of a gricfof 
discrete source values, wherein at least one target value situated among the source values is 
5 determined by means of interpolation. 

Such a method can be used for data of diverse nature, and particularly to expand, compress 
or decompress computer data with a non-critical bit value, such as for instance sound and 
animation files. The starting point is always a one- or multi-dimensional grid of discrete 

10 source values among which concrete target values have to be predicted by means of 

interpolation. The source values herein comprise numerical values which, in the case of for 
instance an image, indicate the colour intensity of a basic colour present therein and, in the 
case of a sound file, represents the frequency, or pitch, or the intensity thereof In the case of 
other types of data a different information component in the source data will similarly be 

1 5 taken as measure . 

For digital images there exists a number of more or less standardized forms in which the 
data is stored. The so-called bitmap format is for instance a form wherein the data is stored 
in a series of discrete numerical values. This series can be converted comparatively simply 
20 into a two-dimensional matrix of source values which each indicate the colour intensity of 
pixels corresponding therewith. In a monochromatic image there thus results a single 
matrix, while in a colour image such a matrix can be constructed for each of the basic 
colours. 

25 One of the simplest ways of enlarging such an image is to expand the discrete pixels. The 
result hereof is shown in figures 1 and 4, wherein the lower image Tl is roughly a five-fold 
enlargement of the source image S. It can be seen clearly that such a method of enlargement 
detracts significantly from the resolution of the final image. Because each pixel of the 
source image S reappears in the result Tl as a rectangular image area of a uniform colour 

30 (intensity), a clear block structure is created which does not correspond with the original 
source image. 
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A better result is obtained by determining target values among the source values in the 
output data by interpolation. In the case of an image the most probable colour value of an 
intermediate target value is herein determined on the basis of the colour values of a group of 
surrounding pixels. The starting point here is preferably a weighted average of the colour 

5 values of the surrounding pixels, wherein a heavier weighting is normally assigned to pixels 
located closer than to pixels located further away. The result of such a method, which is 
applied by many modern graphic computer programs in a form which may differ somewhat 
from each other, is shown in figures 2 and 5. The point of departure is the same source 
image S and the same enlargement factor as in the other figures to obtain the end result T2. 

10 A block structure as in a simple pixel expansion is thus avoided to a significant extent in 

that pixels of the source material S are as it were smeared in the end result T2. While for an 
image such as that in figure 5 this does produce a reasonable result, the image of figure 2 
still loses a lot of its definition because the edges of the image become relatively blurred. 
Such blurred edges are relatively unavoidable assuming an interpolation technique since 

15 interpolation between black and white as in the present image will always produce a certain 
grey value at the edges. The larger the enlargement chosen, the more pronounced this effect 
becomes. 

The present invention has for its object to provide a method of the type stated in the 
20 preamble which to at least a significant extent prevents such a blurring of the image, so that 
the quality of the image is better preserved. 

In order to achieve the stated objective a method of the type stated in the preamble has the 
feature according to the invention that in a region round the target value for detennining a 

25 local minimum and a local maximum is determined of the surrounding source values, that in 
a region round the target value to be determined a measure for the dynamics of the 
surrounding source values is determined and that the target value calculated by interpolation 
is adjusted in the direction of either the local maximum or the local minimum on the basis 
of the calculated measure for the dynamics of the surrounding source values. Not only is an 

30 intermediate value herein determined by interpolation, an analysis is also carried out in the 
source data in order to define the dynamics or contrast thereof in at least the region where a 
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new target value has to be calculated. On the basis of the thus determined dynamics the 
interpolated value is adjusted more or less to the value of a local minimum or maximum. In 
the case of for instance an already relatively blurred source image, the thus adjusted target 
value will not differ from the interpolated value, or hardly so, but in the case of a relatively 
5 dynamic image with sharp local differences in the source values this adjustment will be 
capable of shifting the interpolated value considerably toward the local maximum or 
minimum so that the contrast in the source data is retained to a significant degree. The result 
T3 of an embodiment of the method according to the invention is illustrated in figures 3 and 
6, utilizing the same source image S and expansion factor as in the other figures, and clearly 
10 shows the improvement compared to the existing expansion techniques. 

In a particular embodiment the method according to the invention has the feature that a 
source value which lies closest to the target value for determining in the grid is taken as 
source of a region extending over a finite number of mutually adjacent source values and 

1 5 that the local maximum and the local minimum are determined in this region. In order to 

determine the local minimum and maximum for each target value which has to be calculated 
a uniform framework is herein always used as starting point with just as large a number of 
source values as reference. The same applies for the determining of the local dynamics of 
the source values in a further particular embodiment of the method according to the 

20 invention, characterized in that the measure for the dynamics of the source values is 

determined in a second region extending over a finite number of mutually adjacent source 
values, which second region is optionally of the same size as the first region in which the 
local maximum and minimum are determined. 

25 In the context of the present invention local dynamics in the source values are understood to 
mean the degree in which the source values in the relevant region differ from each other and 
also the steepness with which these differences are present. A further particular embodiment 
of the method according to the invention comprises in this respect an algorithm which 
assigns a numerical value hereto and is characterized for this purpose in that the dynamics 

30 are derived from a normalized difference between a source value and an average of all 

source values in the second region. For the average of all source values in the second region 
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a weighted average is herein preferably taken which assigns a heavier weighting to source 
values located closer in the grid than to source values located further away and which 
particularly utilizes a non-linear density distribution for the purpose of determining the 
weighting factors and more particularly a Gaussian distribution, at least an exponential 
5 density distribution. 

In order to determine the different factors which play a part in predicting the target value 
within the scope of the invention, the starting point is always a number of source values in a 
finite first and second region as frame of reference. Relatively far-removed source values 
can be ignored because they make no or hardly any contribution towards the precision with 
which the target value is calculated and can even detract from this precision in that there is a 
lack of any relationship with the target value for predicting- The calculation time of the 
method moreover increases when more source values are taken as frame of reference. 
Conversely, only one or two reference values will in many cases be too few to enable 
prediction of the target value in reliable manner. A particular embodiment of the method 
according to the invention has been found effective in practice wherein the first and second 
region both extend over nine source values in the grid. 

The adjustment of the interpolated value to the local maximum or minimum within the 
20 scope of the present invention can be performed per se in diverse ways. A further particular 
embodiment of the method according to the invention has the feature however that the final 
target value is a weighted average of the target value determined on the basis of 
interpolation and the local maximum and minimum, wherein a weighting factor is employed 
which depends on average local dynamics of the source values located around the target 
25 value and the relative location of the target value determined on the basis of interpolation 
relative to the local maximum and minimum. 

In order to determine the local minima, maxima and dynamics it is per se possible to carry 
out a separate calculation for each target value to be determined. A preferred embodiment of 
30 the method according to the invention has the feature however that for all source values an 
associated local minimum and maximum and an associated dynamic value is determined 
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beforehand so as to be read for the purpose of determining the target value. A calculation is 
herein performed once for all source values together in order to determine said values, so 
that these are then immediately available. This saves considerable calculation time since 
calculations which would otherwise be performed repeatedly now take up calculation time 
5 once-only. 

The invention is not only suitable for data with only a single information component but 
also for data in the form of source values with separate numerical values for separate 
information components. In accordance with a further embodiment of the method according 
10 to the invention, for each information component a target value located among the source 
values is herein determined individually for this information component. 

The invention will now be further elucidated with reference to an embodiment and an 
associated drawing. In the drawing: 
1 5 figure 1 shows the enlargement of a first image based on an existing technique 

making using of pixel expansion; 
figure 2 shows an enlargement of the first image based on a second existing technique 

making using solely of pixel interpolation; 
figure 3 shows an enlargement of the first image based on an embodiment of the 

20 method according to the present invention; 

figure 4 shows the enlargement of a second image based on an existing technique 

making using of pixel expansion; 
figure 5 shows an enlargement of the second image based on a second existing 

technique making using solely of pixel interpolation; 
25 figure 6 shows an enlargement of the second image based on the embodiment of the 

method according to the present invention used in figure 3. 

Although the invention lends itself in principle to any type of data wherein the exact bit 
value is not critical, it is best elucidated within the scope of the present invention on the 
30 basis of an example with graphic data in the form of an image S. The data on the basis of 
which the image can be constructed is possibly already in the form of a grid or matrix 
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corresponding with the pixels of the image, with for each pixel a numerical source value 
indicating a colour value or, with a view to the method according to the invention, is 
brought into such a form. In the case of a monochromatic image the colour value amounts 
only to an intensity, for instance a grey value; in the case of a colour image this will be a set 
5 of for instance three values, one for each basic colour. A much used format in which images 
are digitally stored is the so-called bitmap format which, for each pixel of a colour image, 
comprises three channels of 8 bits each, and therefore has 256 values for each basic colour. 
This format can be employed relatively directly by the method according to the invention, 
other formats possibly require a conversion similar to that with which such a format is 
10 transferred to the (image) memory of a computer. 

Assuming a black/white image S, a part of the thus assembled matrix is shown by the grid 
of figure 7, wherein the points in the grid represent the pixels of the source image S with for 
each grid point Sy a discrete value of Iy from 0 to 255 from the matrix to thereby designate 
15 the grey value of the relevant pixel. In order to enlarge such an image while retaining 

resolution, intermediate points have to be created, of which the value of the grey value has 
to be calculated. One of these points T is shown in the figure. 

A first estimate of the grey value I t of the intermediate pixel T can be obtained by 
20 interpolation on the basis of the values of the surrounding pixels. The point of departure 
herefor in this example is a weighted average of the values of the surrounding pixels, 
wherein the weighting assigned to each pixel is highly dependent on its distance from T. A 
more particular point of departure herein for the weighting factors is a normalized Gaussian 
distribution, of which pixel T forms the origin. Owing to the sharp fall in the weighting . 
25 factors according to this curve, it is possible to suffice for the interpolation with a finite 

region A of for instance 4x4 pixels around T, which pixels are hatched in the figure. On the 
basis of this interpolation there results an interpolated grey value P t for the pixel T to be 
calculated in accordance with the following relation: 



30 According to the invention a local maximum and a local minimum around the pixel T to be 
calculated is also determined in addition to the interpolated value P r In the embodiment 
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described here the starting point for this purpose is a reference region B of 4x4 pixels S.,_. 
i,S 0>0 S 2 2 , an origin S 0 0 of which is formed by the pixel in the original image lying just in 
front of the pixel for calculating. This region is shown hatched in the figure. A local 
maximum I mfly and a local minimum I mm of the grey values of the original pixels are 
determined within this region. 

In addition, the local dynamics in the grey values are determined according to the invention 
within a reference region around the pixel T for calculating. The starting point herefor in this 
embodiment is the same reference region B as that in which the local minimum and 
maximum have been determined. These dynamics provide a measure for the hardness or 
contrast of the image and are derived in this embodiment from the normalized average 
difference between the grey values of the pixels Sy in reference region B and a weighted 
average of the pixels in a reference region C of 5x5 pixels around the pixel Sy for which the 
local dynamics has to be calculated. This region C therefore differs for each pixel S y and is 
indicated by way of example in the figure for Sqq. This is expressed in the following 
formula: 

(/>)€ =£ V ' || U gem(ij) 



'J 



"here J gem{iJ) = E tMl-c nith^^G^ = 1 



The starting point for the weighted average here are weighting factors G'^ on the basis 
of a normalized Gaussian distribution which attributes a greater weighting to pixels close to 
the pixel T for calculating than to further removed pixels. For the resulting factor Dy there 
applies: O^Dy^ 1 . A value 0^=1 indicates that the source image in the relevant region is very 
hard, i.e. has very sharp transitions in intensity, while a value D^pO indicates precisely the 
opposite, i.e. no differences in intensity occur within the region. The local dynamics D t 
associated with the pixel T to be calculated are derived from the weighted average of the 
values D y from the reference region B around the pixel T to be calculated. 
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The grey value of the pixel T to be calculated is established by adjusting the interpolated 
value P t to the local maximum or minimum on the basis of the thus found hardness value D ( . 
Taken for this purpose in this embodiment is a weighted average of P t and the local 
maximum and minimum, wherein a weighting factor is applied which depends on the 
5 difference from P t and is proportional to the hardness value D t . This is expressed for 
instance by the following formula applied in this example: 

/, = P,.{\ - D n + D,.[I min + (/„„ - / min )-sin„( ^"_ 7 7" )] 

''max ''min 

wheres'm 0 (t) = ! A + ^A.sinfxft-'A) and sm n (t) = sin 0 .sin,^ (t) with n * 5 

If in the relevant reference region B the source image is completely level and has no 
differences in intensity, which is reflected by a value D=0, no shifting or adjustment of the 

10 interpolated value P t takes place. In the case of a comparatively hard image with great 

differences in intensity, wherein D will move in the direction of a value D=l, the adjustment 
will in contrast be maximal. The final value I, will in that case be shifted to the local 
minimum or maximum, depending on which value lies closest to the interpolated value P t: , 
thereby resulting in more definition. The local dynamics of the original image S are thus 

1 5 taken into account in the final grey value I t which is calculated for the pixel T to be added. 

The foregoing algorithm is applied for all pixels which have to be calculated among or 
adjacently of the original pixels in order to create the desired enlargement while retaining 
resolution. In practice this means that for all pixels of the original source image the 

20 associated local maxima, minima and average hardness values are calculated beforehand so 
that these are then immediately available for the above stated calculations. The stated 
algorithm allows of simple translation into computer software with which a suitable 
computer can be loaded to perform the calculations fully automatically. In a user interface 
adapted thereto setting options can be offered for fme-runing of the size of the reference 

25 regions, the relations for the weighting factors and other parameters. 

The result on the basis of this embodiment of the method according to the invention is 
shown in figures 3 and 6 for respectively a relatively hard source image in the form of a 
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sharply defined character and a more variegated source image corresponding with a 
photographic image. It can clearly be seen that the definition of the character in figure 3 has 
been retained despite the great enlargement thereof, while the variegation in the image of 
figure 6 is retained with the method according to the invention and is even better preserved 
5 than on the basis of solely an interpolation method as shown in figure 5. 

Although the invention has been further elucidated above with reference to only two 
embodiments, it will be apparent that the invention is by no means limited to the given 
embodiments. The invention can thus also be applied to images with more colours, wherein 

10 the above stated algorithm is applied separately to each basic colour. The invention can 

advantageously also be applied for the interpolation of sound data and other data with non- 
critical bit values while retaining dynamics. The above stated weighting factors and frames 
of reference, although very effective, have been given solely as examples. Depending on the 
concrete situation, a fine-tuning thereof can take place to further enhance the quality of the 

15 final result. Larger, smaller or identical reference regions can for instance thus be applied to 
calculate the various above stated factors and, in order to determine the final grey value of 
the pixel for calculating, a different algorithm can also be chosen which takes account of the 
local dynamics derived in this pixel and the local maximum and minimum. 

20 The invention generally provides a method of interpolating data while retaining dynamics, 
which implies that the character and the definition of the source data is retained to a 
significant degree in the result. 
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Claims 

1. Method for processing data in the form of a grid of discrete source values, wherein at 
least one target value situated among the source values is determined by means of 

5 interpolation, characterized in that in a region round the target value for detennining a 

local minimum and a local maximum is determined of the surrounding source values, that in 
a region round the target value to be determined a measure for the dynamics of the 
surrounding source values is determined and that the target value calculated by interpolation 
is adjusted in the direction of either the local maximum or the local minimum on the basis of 
10 the calculated measure for the dynamics of the surrounding source values. 

2. Method as claimed in claim 1, characterized in that use is made of interpolation on 
the basis of a non-linear density distribution which assigns a heavier weighting to source 
values located closer in the grid than to source values located further away, in particular a 

1 5 Gaussian distribution, at least an exponential density distribution. 

3. Method as claimed in claim 1 or 2, characterized in that a source value which lies 
closest to the target value for determining in the grid is taken as source of a region extending 
over a finite number of mutually adjacent source values and that the local maximum and the 

20 local minimum are determined in this region. 

4. Method as claimed in claim 3, characterized in that the measure for the dynamics 
of the source values is determined in a second region extending over a finite number of 
mutually adjacent source values, which second region is optionally of the same size as the 

25 first region in which the local maximum and minimum are determined. 

5. Method as claimed in claim 4, characterized in that the dynamics are derived from 
a normalized difference between a source value and an average of all source values in the 
second region. 



30 
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6. 



Method as claimed in claim 5, characterized in that for the average of all source 



values in the second region a weighted average is taken which assigns a heavier weighting to 
source values located closer in the grid than to source values located further away and which 
particularly utilizes a non-linear density distribution for the purpose of determining the 
5 weighting factors and more particularly from a Gaussian distribution, at least from an 
exponential density distribution. 

7. Method as claimed in any of the foregoing claims, characterized in that the final 
target value is a weighted average of the target value determined on the basis of 
1 0 interpolation and the local maximum and minimum, wherein a weighting factor is employed 
which depends on average local dynamics of the source values located around the target 
value to be determined and the relative location of the target value determined on the basis 
of interpolation relative to the local maximum and minimum. 

15 8. Method as claimed in any of the foregoing claims, characterized in that for all 

source values an associated local minimum and maximum and an associated dynamic value 
are determined beforehand so as to be read for the purpose of determining the target value. 

9. Method as claimed in any of the foregoing claims, characterized in that the source 
20 values comprise separate numerical values for separate information components in the 

source data and that for each information component individually a target value located 
among the source values is determined for this information component. 

10. Method as claimed in claim 9, characterized in that the data comprises the pixels 
25 of an image with a separate discrete source value for each basic colour therein and that for 

each basic colour individually a target value located among the source values is determined. 

11. A computer loaded with software for performing the method as claimed in one or 
more of the foregoing claims. 



30 



WO 00/73993 ^^^^ ^^^^ PCT/NLO 0/00359 

12 

12. Software for performing the method as claimed in one or more of the claims 1 - 1 0 in 
a computer. 
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Werkwijze voor het verwerken van data. 



De onderhavige uitvinding heeft betrekking op een werkwijze voor het verwerken van 
data in de vorm van een rooster van discrete bronwaarden, waarbij door middel van 
5 interpolate ten minste een tussen de bronwaarden gelegen doelwaarde wordt bepaald. 

Een dergelijke werkwijze is bruikbaar voor data van uiteenlopende aard en met name 
voor het vergroten, verdichten of de-comprimeren van computerdata met een niet- 
kritische bitwaarde, zoals bijvoorbeeld geluids- en animatiebestanden. Het uitgangspunt 

10 is steeds een een- of meerdimensionaal rooster van discrete bronwaarden waartussen 
door middel van interpolate concrete doelwaarden dienen te vvorden voorspeld. De 
bronwaarden omvatten daarbij getalswaarden die in het geval van bijvoorbeeld een 
afbeelding de kleurintensiteit van een daarin aanwezige basiskleur aangeeft en in het 
geval van een geluidsbestand de frequence, oftewel toonhoogte, dan wel de intensiteit 

15 daarvan representeert. In het geval van andersoortige data zal overeenkomstig een 
andere informatiecomponent in de uitgangsdata als maat worden genomen. 

Voor digitate afbeeldingen bestaat een aantal min of meer gestandaardiseerde vormen 
waarin de data worden opgeslagen. Het zogenaamde bitmap formaat is bijvoorbeeld een 

20 vorm waarbij de data in een reeks van discrete getalswaarden liggen opgeslagen. Deze 
reeks kan betrekkelijk eenvoudig worden omgezet in een tweedimensionale matrix van 
bronwaarden die ieder de kleurintensiteit van daarmee corresponderende beeldpunten 
aangeven. Bij een monochromatische afbeelding ontstaat aldus een enkele matrix, 
terwijl bij een kleurenafbeelding een dergelijke matrix kan worden geconstrueerd voor 

25 ieder van de basiskleuren. 

Een van de meest eenvoudige wijzen om een dergelijke afbeelding te vergroten is het 
vergroten van de discrete beeldpunten. Het resultaat daarvan is in figuur 1 en 4 
aangeven, waarbij de onderste afbeelding Tl een circa vijfvoudige uitvergroting is van 
30 de uitgangsafbeelding S. Duidelijk is te zien dat een dergelijke vvijze van uitvergroting 
in belangrijke mate afbreuk doet aan de resolutie van de uiteindelijke afbeelding. 
Doordat ieder beeldpunt van de uitgangsafbeelding S in het resultaat Tl als een 
rechthoekig beeldvlak van een uniforme kleur(intensiteit) terugkomt, ontstaat een 
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duidelijke blokstructuur die niet overeenstemt met de oorspronkelijke 
uitgangsafbeelding. 

Een beter resultaat wordt verkregen door tussen de bronwaarden in de uitgangsdata 

5 doelwaarden door interpolatie vast te stellen. In het geval van een afbeelding wordt 

daarbij op basis van de kleurwaarden van een groep omringende beeldpunten de meest 
waarschijnlijke kleurwaarde van een tussengelegen doelwaarde bepaald. Bij voorkeur 
wordt daarbij uitgegaan van een gewogen gemiddelde van de kleurwaarden van de 
omringende beeldpunten, waarbij gewoonlijk aan meer nabijgelegen beeldpunten een 

10 zwaarder gewicht wordt toegekend dan aan verder weg gelegen beeldpunten. Het 
resultaat van een dergelijke methode die door veel hedendaagse grafische 
computerprogramma's in onderling al of niet ietwat verschillende vorm wordt toegepast 
is in figuur 2 en 5 weergegeven. Uitgegaan is van dezelfde uitgangsafbeelding S en 
eenzelfde vergrotingsfactor als in de andere figuren ter verkrijging van het eindresultaat 

15 T2. Een blokstructuur zoals bij een simpele beeldpuntvergroting wordt aldus in 

belangrijke mate vermeden doordat beeldpunten van het bronmateriaal S als het ware 
zijn uitgesmeerd in het eindresultaat Tl. Voor een afbeelding zoals die van figuur 5 
levert dit weliswaar een redelijk resultaat, doch de afbeelding van figuur 2 verliest veel 
van diens scherpte doordat de randen van de afbeelding relatief vaag raken. Dergelijke 

20 vage randen zijn betrekkelijk onvermijdelijk uitgaande van een interpolatietechniek 

omdat interpolatie tussen zwart en wit zoals in de onderhavige afbeelding aan de randen 

i 

steeds een zekere grijswaarde zal opleveren. Naarmate een grotere uitvergroting wordt 
gekozen wordt dit effect sterker. 

25 Met de onderhavige uitvinding wordt beoogd in een vverkwijze van de in de aanhef 
genoemde soort te voorzien die een dergelijke vervaging van de afbeelding althans in 
belangrijke mate tegengaat, opdat de kwaliteit van de afbeelding beter wordt behouden. 



30 



Om het beoogde doel te bereiken heeft een werkwijze van de in de aanhef genoemde 
soort volgens de uitvinding als kenmerk dat in een gebied om de te bepalen doelwaarde 
een lokaal minimum en een lokaal maximum wordt bepaald van de omringende 
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bronvvaarden, dat in een gebied om de te bepalen doelwaarde een maat voor de 
dynamiek van de omringende bronwaarden wordt bepaald en dat de door interpolate 
berekende doelwaarde op basis van de berekende maat voor de dynamiek van de 
omliggende bronwaarden in de richting van hetzij het lokale maximum, hetzij het lokale 



Het resultaat T3 van een uitvoeringsvorm van een de werkwijze volgens de uitvinding is 
in figuur 3 en 6 weergegeven, uitgaande van dezelfde bronafbeelding S en 
vergrotingsfactor als bij de andere figuren en geeft duidelijk de verbetering aan ten 
opzichte van de bestaande vergrotingstechnieken. 

20 

In een bijzondere uitvoeringsvorm heeft de werkwijze volgens de uitvinding als 
kenmerk dat een bronwaarde die in het rooster het meest nabij de te bepalen doelwaarde 
ligt als oorsprong wordt genomen van een gebied dat zich over een eindig aantal 
onderling naburige bronwaarden uitstrekt en dat in dit gebied het lokale maximum en 

25 het lokale minimum wordt bepaald. Hierbij wordt ter bepaling van het lokale minimum 
en maximum voor iedere doelwaarde die dient te worden berekend steeds van een 
uniform kader uitgegaan met een evengroot aantal bronwaarden als referentie. Hetzelfde 
geldt voor de bepaling van de lokale dynamiek van de bronwaarden in een verdere 
bijzondere uitvoeringsvorm van de werkwijze volgens de uitvinding, gekenmerkt 

30 doordat in een tweede gebied dat zich over een eindig aantal onderling naburige 

bronwaarden uitstrekt de maat voor de dynamiek van de bronwaarden wordt bepaald, 



10 



15 



5 



minimum wordt bijgesteld. Hierbij wordt niet alleen door interpolatie een 
tussenliggende waarde bepaald, maar tevens een analyse uitgevoerd in de uitgangsdata 
om de dynamiek of contrast daarvan te bepalen in althans het gebied waar een nieuwe 
doelwaarde dient te worden berekend. Op basis van de aldus bepaalde dynamiek wordt 
de geinterpoleerde waarde meer of minder naar de waarde van een lokaal minimum of 
maximum bijgesteld. In het geval van bijvoorbeeld een reeds betrekkelijk vage 
bronafbeelding zal de aldus bijgestelde doelwaarde niet of nauwelijks van de 
geinterpoleerde waarde afwijken, doch in het geval van een betrekkelijk dynamische 
afbeelding met lokaal scherpe verschillen in de bronwaarden, zal deze bijstelling de 
geinterpoleerde waarde sterk naar het lokale maximum of minimum kunnen 
verschuiven zodat het contrast in de uitgangsdata in belangrijke mate behouden blijft. 
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welk tweede gebied al of niet van eenzelfde grootte is als het eerste gebied waarin het 
lokale maximum en minimum worden bepaald. 

In het kader van de onderhavige uitvinding wordt onder lokale dynamiek in de 
bronvvaarden verstaan de mate waarin de bronwaarden in het betreffende gebied van 
elkaar verschillen en ook de steilheid waarmee deze verschillen aanwezig zijn. Een 
verdere bijzondere uitvoeringsvorm van de werkwijze volgens de uitvinding omvat in 
dit verband een algoritme dat hieraan een getalswaarde toekent en is daartoe gekenmerkt 
doordat de dynamiek wordt gevonden uit een genormeerd verschil van een bronwaarde 
met een gemiddelde van alle bronwaarden in het tweede gebied. Bij voorkeur wordt 
daarbij voor het gemiddelde van alle bronwaarden in het tweede gebied een gewogen 
gemiddelde genomen dat aan meer nabij gelegen bronwaarden in het rooster een 
zwaarder gewicht toekent dan aan verder gelegen bronwaarden en dat in het bijzonder 
uitgaat van een niet lineaire dichtheidsverdeling ter bepaling van de gewichtsfactoren en 
meer in het bijzonder van een Gauss-verdeling, althans een exponentiele 
dichtheidsverdeling. 

Voor het bepalen van de verschillende factoren die binnen het kader van de uitvinding 
een rol spelen voor het voorspellen van de doehvaarde wordt steeds uitgegaan van een 
aantal bronwaarden in een eindig eerste en tweede gebied als referentiekader. Relatief 
ver verwijderd gelegen bronwaarden kunnen buiten beschouwing blijven omdat zij niet 
of nauwelijks bijdragen aan de precisie waarmee de doehvaarde wordt berekend en zelfs 
afbreuk kunnen doen aan die precisie doordat het verband met de te voorspellen 
doehvaarde ontbreekt. Bovendien neemt de rekentijd van de werkwijze toe naarmate 
meer bronwaarden als referentiekader worden genomen. Anderzijds zullen slechts een 
of twee referentiewaarden in veel gevallen te weinig zijn om de doehvaarde op een 
betrouwbare wijze te kunnen voorspellen. In de praktijk is een bijzondere 
uitvoeringsvorm van de werkwijze volgens de uitvinding probaat gebleken, waarbij het 
eerste en tweede gebied zich beiden over negen bronwaarden in het rooster uitstrekken. 
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De bijstelling van de gemterpoleerde vvaarde naar het lokale maximum of minimum 
binnen het kader van de onderhavige uitvinding kan op zichzelf op verscheidene wijzen 
worden uitgevoerd. Een verdere bijzondere uitvoeringsvorm van de werkwijze volgens 
de uitvinding heeft evenwel als kenmerk dat de uiteindelijke doelwaarde een gevvogen 
5 gemiddelde is van de op basis van interpolate bepaalde doelwaarde en het lokale 

maximum en minimum, vvaarbij een gewichtsf actor wordt gehanteerd die afhangt van 
een gemiddelde lokale dynamiek van de rondom de te bepalen doelwaarde gelegen 
bronvvaarden en de relatieve ligging van de op basis van interpolate bepaalde 
doelwaarde ten opzichte van het lokale maximum en minimum. 

10 

Op zichzelf kan ter bepaling van de lokale minima, maxima en dynamiek telkens voor 
iedere te bepalen doelwaarde afzonderlijk een berekening worden uitgevoerd. Een 
voorkeursuitvoeringsvorm van de werkwijze volgens de uitvinding heeft evenwel als 
kenmerk dat voor alle bronwaarden een bijbehorend lokaal minimum en maximum, . 

15 alsmede een bijbehorende dynamiekwaarde vooraf wordt bepaald om ter bepaling van 
de doelwaarde te worden uitgelezen. Hierbij wordt eenmaal voor alle bronwaarden 
tegelijk een berekening uitgevoerd ter bepaling van de genoemde waarden, zodat deze 
vervolgens rechtstreeks beschikbaar zijn. Dit bespaart aanzienlijk in rekentijd doordat 
berekeningen die anders herhaaldelijk zouden worden uitgevoerd nu slechts eenmalig 

20 rekentijd in beslag nemen. 

i 

De uitvinding leent zich niet alleen voor data met slechts een enkele 
informatiecomponent maar ook voor data in de vorm van bronwaarden met 
afzonderlijke getalswaarden voor afzonderlijke informatiecomponenten. Conform een 
25 verdere uitvoeringsvorm van de werkwijze volgens de uitvinding wordt daarbij voor 
iedere informatiecomponent afzonderlijk een tussen de bronwaarden gelegen 
doelwaarde voor die informatiecomponent bepaald. 

De uitvinding zal thans nader worden toegelicht aan de hand van een 
30 uitvoeringsvoorbeeld en een bijbehorende tekening. In de tekening toont: 
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figuur 1 de uitvergroting van een eerste afbeelding op basis van een bestaande 

techniek gebruik makend van beeldpuntvergroting; 
figuur 2 een uitvergroting van de eerste afbeelding op basis van een tweede 

bestaande techniek, gebruik makend van louter beeldpuntinterpolatie; 
5 figuur 3 een uitvergroting van de eerste afbeelding op basis van een 

uitvoeringsvorm van de werkwijze volgens de onderhavige uitvinding; 
figuur 4 de uitvergroting van een tweede afbeelding op basis van een bestaande 

techniek gebruik makend van beeldpuntvergroting; 
figuur 5 een uitvergroting van de tweede afbeelding op basis van een tweede 

10 bestaande techniek, gebruik makend van louter beeldpuntinterpolatie; 

figuur 6 een uitvergroting van de tweede afbeelding op basis van de 

uitvoeringsvorm van de werkwijze volgens de onderhavige uitvinding 

gebruikt in figuur 3. 



15 Hoewel de uitvinding zich in beginsel voor iedere soort van data, waarbij de exacte 

bitwaarde niet kritisch is, laat zij zich binnen het bestek van de onderhavige aanvrage 
het best toelichten aan de hand van een voorbeeld met grafische data in de vorm van een 
afbeelding S. De data op basis waarvan de afbeelding kan worden geconstrueerd, 
bevindt zich mogelijk reeds in een vorm van een rooster of matrix corresponderend met 

20 de beeldpunten van de afbeelding met voor ieder beeldpunt een getalsmatige 

bronwaarde die een kleurvvaarde aangeeft of wordt met het oog op de werkwijze volgens 
de uitvinding in een dergelijke vorm gebracht. In het geval van een monochromatische 
afbeelding behelst de kleurwaarde louter een intensiteit, bijvoorbeeld een grijswaarde; 
bij een kleurenafbeelding zal dit een verzameling van bijvoorbeeld drie waarden zijn, 

25 voor ieder basiskleur een. Een veel gebruikt formaat waarin afbeeldingen digitaal 

worden opgeslagen, is het zogenaamde bitmap formaat dat voor ieder beeldpunt van een 
kleurenafbeelding drie kanalen omvat van ieder 8 bit en derhalve voor ieder basiskleur 
256 waarden kent. Dit formaat kan betrekkelijk rechtstreeks door de werkwijze volgens 
de uitvinding worden aangewend, andere formaten vergen mogelijk een conversie gelijk 

30 aan die waarmee een dergelijk formaat in het (beeld)geheugen van een computer wordt 
overgebracht. 
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Uitgaande van een zwart/wit afbeelding S wordt een deel van de aldus samengestelde 
matrix weergegeven door het rooster van figuur 7, waarbij de punten in het rooster de 
beeldpunten van de uitgangs afbeelding S representeren met voor ieder roosterpunt 
een discrete waarde Iy van 0 tot en met 255 uit de matrix om daarmee de grijswaarde 
5 van het betreffende beeldpunt aan te.geven. Om een dergelijke afbeelding met behoud 
van resolutie te vergroten, dienen tussenliggende punten te worden gecreeerd waarvan 
de waarde van de grijswaarde dient te worden berekend. Een van deze punten T is in de 
figuur aangegeven. 

10 Een eerste schatting van de grijswaarde I t van het tussengelegen beeldpunt T kan door 

interpolate worden verkregen op basis van de waarden van de omringende beeldpunten. 
Hiervoor wordt in dit voorbeeld uitgegaan van een gewogen gemiddelde van de 
waarden van de omringende beeldpunten, waarbij het gewicht dat aan ieder beeldpunt 
wordt toegekend in sterke mate afhangt van diens afstand tot T. Meer in het bijzonder 

15 wordt daarbij voor de gewichtsfactoren uitgegaan van een genormeerde Gaussverdeling, 
waarvan beeldpunt T de oorsprong vormt. Door de sterke daling van de 
gewichtsfactoren volgens deze curve kan voor de interpolate worden volstaan met een 
eindig gebied A van bijvoorbeeld 4x4 beeldpunten rondom T, welke punten in de figuur 
zijn gearceerd. Op basis van deze interpolatie volgt een geinterpoleerde grijswaarde P t 

20 voor het te berekenen beeldpunt T conform de navolgende betrekking: 

Volgens de uitvinding wordt naast de geinterpoleerde waarde P t tevens een lokaal 
maximum en een lokaal minimum rondom het te berekenen beeldpunt T bepaald. In het 
hier beschreven voorbeeld wordt daartoe uitgegaan van een referentiegebied B van 4x4 
beeldpunten S. L . 1( S 0i0 ..S 2i2 waarvan een oorsprong S 00 wordt gevorrnd door het beeldpunt 
25 dat in de oorspronkelijke afbeelding juist voor het te berekenen beeldpunt ligt. Dit 
gebied is in de figuur gearceerd weergegeven. Binnen dit gebied wordt een lokaal 
maximum l m:tx en een lokaal minimum 1^ van de grijs waarden van de oorspronkelijke 
beeldpunten bepaald. 
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Voorts wordt volgens de uitvinding binnen een referentiegebied rondom het te 
berekenen beeldpunt T de lokale dynamiek in de grijsvvaarden bepaald. Hiertoe wordt in 
dit voorbeeld uitgegaan van hetzelfde referentiegebied B als dat waarin het lokale 
minimum en maximum werden bepaald. Deze dynamiek geeft een maat voor de 
hardheid oftewel het contrast van de afbeelding en wordt in dit voorbeeld gevonden uit 
het genormeerde gemiddelde verschil tussen de grijswaarden van de beeldpunten Sij in 
het referentiegebied B en een gewogen gemiddelde van de beeldpunten in een 
referentiegebied C van 5x5 beeldpunten rondom het punt S i} waarvoor de lokale 
dynamiek dient te vvorden berekend. Dit gebied C verschilt dus voor ieder punt Sy en is 
in de figuur bij wijze van voorbeeld voor Sqq aangegeven. Een en ander komt tot 
uitdrukking in de navolgende formule: 



•J 



Hierbij wordt voor het gewogen gemiddelde I gem uitgegaan van gewichtsfactoren G' pq op 
basis van een genormeerde Gaussverdeling die aan beeldpunten nabij het te berekenen 
beeldpunt T een groter gewicht hecht dan aan verder verwijderde beeldpunten. Voor de 
resulterende factor geldt: 0<Dij< 1. Een waarde D^I geeft aan dat de bronafbeelding 
in het betreffende gebied zeer hard is, dat wil zeggen zeer scherpe overgangen in 
intensiteit kent, terwijl een waarde D^O juist het tegenovergestelde aangeeft, namelijk 
dat geen intensiteitsverschillen binnen het gebied voorkomen. De met het te berekenen 
beeldpunt T geassocieerde lokale dynamiek D, wordt gevonden uit het gewogen 
gemiddelde van de waarden Djj uit het referentiegebied B rondom het te berekenen punt 
T. 



De grijs waarde van het te berekenen beeldpunt T wordt vastgesteld door op basis van de 
aldus gevonden hardheids waarde D t de geinterpoleerde waarde P t naar het lokale 
maximum of minimum bij te stellen. Hiertoe wordt in dit voorbeeld een gewogen 
gemiddelde genomen van P t en het lokale maximum en minimum, waarbij een 
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weegfactor wordt toegepast die afhangt van het verschil met Pt en evenredig is met de 
hardheidswaard D t . Een en ander wordt bijvoorbeeld tot uitdrukking gebracht door de 
navolgende, in dit voorbeeld toegepaste, formule: 



/, = P,.(l - D) + D,[/, nin+ (/_- / min ).sin„( ' _ )] 

max min 

waarbij sin 0 (f) = Vi + V2.s\n{it{t-Vi)) en sin /j (r)=sin 0 .sin ^jCf) met n~5 

5 Indien in het betreffende referentiegebied B de bronafbeelding volstrekt vlak is en geen 
intensiteitsverschillen kent, vvat wordt weerspiegeld door een waarde D=0, vindt geen 
verschuiving of bijstelling van de geinterpoleerde waarde P t plaats, Bij een betrekkelijk 
harde afbeelding met grote intensiteitsverschillen, waarbij D in de richting van een 
waarde D=l zal gaan, zal de bijstelling daarentegen juist maximaal zijn. De 

10 uiteindelijke waarde I t zal in dat geval naar het lokale minimum of maximum worden 
getrokken, afhankelijk welke waarde het dichts bij de geinterpoleerde waarde P t ligt, 
waardoor meer scherpte ontstaat. Aldus wordt de lokale dynamiek van de 
oorspronkelijke afbeelding S verdisconteerd in de uiteindelijke grijswaarde I, die voor 
het toe te voegen punt T wordt berekend. 

15 

Het voorgaande algoritme wordt toegepast voor alle beeldpunten die tussen of naast de 
oorspronkelijke beeldpunten dienen te worden berekend om de gewenste uitvergroting 
met behoud van resolutie te creeren. In de praktijk betekent dit dat vooraf voor alle 
beeldpunten van de oorspronkelijke bronafbeelding de daarbij behorende lokale 

20 maxima, minima en gemiddelde hardheidswaarden worden berekend, opdat die 

vervolgens rechtstreeks beschikbaar zijn voor de hiervoor aangegeven berekeningen. 
Het aangegeven algoritme laat zich eenvoudig vertalen in computerprogrammatuur 
waarmee een geschikte computer kan worden geladen om de berekeningen volledig 
automatisch uit te voeren. In een daarop afgestemde gebruikers interface kunnen 

25 instellingsmogelijkheden worden geboden om de grootte van de referentiegebieden, de 
relaties voor de gewichtfactoren en andere parameters fijn te regelen. 
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Het resultaat op grond van dit uitvoeringsvoorbeeld van de werkwijze volgens de 
uitvinding is in figuur 3 en 6 weergegeven voor respectievelijk een relatief harde 
bronafbeelding in de vorm van een scherp afgebakend karakter en een meer 
geschakeerde bronafbeelding overeenkomend met een fotografische opname. Duidelijk 
5 is te zien dat de scherpte van het karakter in figuur 3 behouden blijft ondanks de sterke 
uitvergroting daarvan, terwijl ook de schakering in de afbeelding van figuur 6 bij de 
werkwijze volgens de uitvinding behouden blijft en zelfs beter wordt gehandhaafd dan 
op grond van louter een interpolatiemethode zoals weergegeven in figuur 5. 

10 Hoewel de uitvinding hiervoor aan de hand van louter een tweetal voorbeelden nader 
werd toegelicht zal het duidelijk zijn dat de uitvinding geenszins tot het gegeven 
voorbeeld is beperkt. Zo kan de uitvinding ook worden toegepast op afbeeldingen met 
meer kleuren, waarbij het hiervoor aangegeven algoritme op ieder basiskleur 
afzonderlijk wordt toegepast. Ook kan de uitvinding met voordeel worden toegepast 

15 voor het met behoud van dynamiek interpoleren van geluidsdata en andere data met 

niet-kritische bitwaarden. De hiervoor aangegeven gewichtsfactoren en referentiekaders 
werden, hoewel zeer probaat, louter als voorbeeld gegeven. Afhankelijk van een 
concrete situatie kan daarvoor een fijnaf stemming plaatsvinden om de kwaliteit van het 
eindresultaat verder te bevorderen. Zo kunnen bijvoorbeeld grotere. kleinere of 

20 identieke referentiegebieden worden toegepast voor het berekenen van de verscheidene 
hiervoor aangegeven factoren en ook kan voor het bepalen van de uiteindelijke 
grijswaarde van het te berekenen punt een ander algoritme worden gekozen dat rekening 
houdt met de gevonden Iokale dynamiek in dit punt en het lokale maximum en 
minimum. 

25 

In het algemeen biedt de uitvinding een werkwijze om data met behoud van dynamiek te 
interpoleren, wat inhoudt dat het karakter en de scherpte van de uitgangsdata in 
belangrijke mate wordt gehandhaafd in het resultaat. 
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Conclusies: 

1 . Werkwijze voor het verwerken van data in de vorm van een rooster van discrete 
bronwaarden, waarbij door middel van interpolatie ten minste een tussen de 

5 bronwaarden gelegen doelwaarde wordt bepaald, met het kenmerk dat in een gebied om 
de te bepalen doelwaarde een lokaal minimum en een lokaal maximum wordt bepaald 
van de omringende bronwaarden, dat in een gebied om de te bepalen doelwaarde een 
maat voor de dynamiek van de omringende bronwaarden wordt bepaald en dat de door 
interpolatie berekende doelwaarde op basis van de berekende maat voor de dynamiek 

10 van de omliggende bronwaarden in de richting van hetzij het lokale maximum, hetzij 
het lokale minimum wordt bijgesteld. 

2. Werkwijze volgens conclusie 1 met het kenmerk dat gebruik wordt gemaakt van 
interpolatie op basis van een niet-lineaire dichtheidsverdeling die aan meer nabij 

15 gelegen bronwaarden in het rooster een zwaarder gewicht toekent dan aan verder 

gelegen bronwaarden, in het bijzonder een Gauss-verdeling, althans een exponentiele 
dichtheidsverdeling. 

3. Werkwijze volgens conclusie 1 of 2 met het kenmerk dat een bronwaarde die in 
20 het rooster het meest nabij de te bepalen doelwaarde ligt als oorsprong wordt genomen 

van een gebied dat zich over een eindig aantal onderling naburige bronwaarden uitstrekt 
en dat in dit gebied het lokale maximum en het lokale minimum wordt bepaald. 

4. Werkwijze volgens conclusie 3 met het kenmerk dat in een tweede gebied dat 
25 zich over een eindig aantal onderling naburige bronwaarden uitstrekt de maat voor de 

dynamiek van de bronwaarden wordt bepaald, welk tweede gebied al of niet van 
eenzelfde grootte is als het eerste gebied waarin het lokale maximum en minimum 
worden bepaald. 
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5. Werkwijze volgens conclusie 4 met het kenmerk dat de dynamiek wordt 
gevonden uit een genormeerd verschil van een bronwaarde met een gemiddelde van alle 
bronwaarden in het tvveede gebied. 

5 6. Werkwijze volgens conclusie 5 met het kenmerk dat voor het gemiddelde van 

alle bronwaarden in het tweede gebied een gewogen gemiddelde wordt genomen dat aan 
meer nabij gelegen bronwaarden in het rooster een zwaarder gewicht toekent dan aan 
verder gelegen bronwaarden en in het bijzonder uitgaat van een niet lineaire 
dichtheidsverdeling ter bepaling van de gewichtsfactoren en meer in het bijzonder van 
10 een Gauss-verdeling, althans van een exponentiele dichtheidsverdeling. 

7. Werkwijze volgens een der voorafgaande conclusies met het kenmerk dat de 
uiteindelijke doelwaarde een gewogen gemiddelde is van de op basis van interpolate 
bepaalde doelwaarde en het lokale maximum en minimum, waarbij een gewichtsfactor 

15 wordt gehanteerd die afhangt van een gemiddelde lokale dynamiek van de rondom de te 
bepalen doelwaarde gelegen bronwaarden en de relatieve ligging van de op basis van 
interpolate bepaalde doelwaarde ten opzichte van het lokale maximum en minimum. 

8. Werkwijze volgens een der voorafgaande conclusies met het kenmerk dat voor 
20 alle bronwaarden een bijbehorend Iokaal minimum en maximum, alsmede een 

bijbehorende dynamiek waarde vooraf wordt bepaald om ter bepaling van de doelwaarde 
te worden uitgelezen. 

9. Werkwijze volgens een der voorafgaande conclusies met het kenmerk dat de 
25 bronwaarden afzonderlijke getalswaarden omvatten voor afzonderlijke 

informatiecomponenten in de uitgangsdata en dat voor iedere informatiecomponent 
afzonderlijk een tussen de bronwaarden gelegen doelwaarde voor die 
informatiecomponent wordt bepaald. 



30 



10. Werkwijze volgens conclusie 9 met het kenmerk dat de data de beeldpunten van 
een afbeelding omvatten met voor iedere basiskleur daarin een afzonderlijke discrete 
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bronwaarde en dat voor iedere basiskleur afzonderlijk een tussen de bronvvaarden 
gelegen doelwaarde wordt bepaald. 

1 1 . Een computer geladen met programmatuur voor het uitvoeren van de werkvvijze 
5 volgens een of meer der voorafgaande conclusies. 

12. Programmatuur voor het uitvoeren van de werkwijze volgens een of meer der 
conclusies 1 tot en met 10 in een computer. 



t 
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Uittreksel: 

Werkvvijze voor het verwerken van data. 

5 Bij een werkwijze voor het verwerken van data in de vorm van een rooster van discrete 
bronwaarden wordt door middel van interpolate ten minste een tussen de bronwaarden 
gelegen doelwaarde bepaald. Volgens de uitvinding wordt tevens rekening gehouden 
met het karakter en de scherpte van de brondata door in een gebied om de te bepalen 
doelwaarde een lokaal minimum en een lokaal maximum van de omringende 

10 bronwaarden en in een gebied om de te bepalen doelwaarde een maat voor de dynamiek 
van de omringende bronwaarden te bepalen. De door interpolate berekende doelwaarde 
wordt op basis van de aldus berekende maat voor de dynamiek van de omliggende 
bronwaarden in de richting van hetzij het lokale maximum, hetzij het lokale minimum 
bijgesteld om uitdrukking te geven aan het karakter en de scherpte van de uitgangsdata. 

15 



Fig. 3 



